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ABSTRACT 

We present simulations of Triton's post-capture orbit that confirm the im- 
portance of Kozai-type oscillations in its orbital elements. In the context of 
the tidal orbital evolution model, these variations require average pericenter dis- 
tances much higher than previously published, and the timescale for the tidal 
orbital evolution of Triton becomes longer than the age of the Solar System. 
Recently-discovered irregular satellites present a new constraint on Triton's or- 
bital history. Our numerical integrations of test particles indicate a timescale 
for Triton's orbital evolution to be less than 10 5 yrs for a reasonable number of 
distant satellites to survive Triton's passage. This timescale is inconsistent with 
the exclusively tidal evolution (time scale of > 10 8 yrs), but consistent with the 
interestion with the debris from satellite-satellite collisions. Any major regular 
satellites will quickly collide among themselves after being perturbed by Triton, 
and the resulting debris disk would eventually be swept up by Triton; given that 
the total mass of the Uranian satellite system is 40% of that of Triton, large scale 
evolution is possible. This scenario could have followed either collisional or the 
recently-discussed three-body-interaction-based capture. 

Subject headings: celestial mechanics — planets and satellites: individual(Neptune, Triton) 



1. Introduction 

Neptune's Triton is the only large planetary satellite to orbit retrograde relative to the 
planet's rotation. McCord (1966) and later McKinnon (1984) suggested that Triton is a 
captured satellite, whose originally eccentric orbit was circularized due to tidal dissipation 
within Triton. Goldreich et al. (1989) postulate that Triton was captured from heliocentric 
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orbit by a collision with a pre-existing satellite, and its initial high-eccentricity orbit then 
evolved due to tidal dissipation alone. They showed that the tidal evolution timescale is 
significantly shorter than the age of the Solar System (a few times 10 s years), even when 
the variations in Triton's eccentricity with 1/2-yr period are accounted for. However, McK- 
innon et al. (1995) noted that the criteria for non- disruptive capture are much stricter than 
Goldreich et al. (1989) calculated. If the original regular satellites of Neptune and Uranus 
were similar, a collision with the largest moons (preferred due to their cross-sections) would 
disrupt Triton, with re-accretion on an orbit inclined to Neptune's equator being impossible. 

McKinnon & Leith (1995) suggested instead that Triton was captured and its orbit was 
evolved by aerodynamic drag in Neptune's primordial protosatellite nebula, and that after 
its orbit ewas circularized Triton's gravity would be sufficient to clear a disk gap and thus 
halt further evolution. Gas drag has been suspected as a capture mechanism for small and 
distant irregular satellites, but capture of Triton would require unprecedented gas densities, 
requiring very close approaches to Neptune. 

Agnor & Hamilton (2004) propose a three-body capture scenario for Triton. They 
suggest pre-capture Triton may have been a member of a binary, whose disruption during a 
Neptune encounter led to Triton's capture and its companion's escape. Their work addresses 
only the capture itself, leaving the physics of post-capture evolution unchanged. 

2. Effect of Solar Perturbations 

Any capture mechanism, be it collision, gas drag or 3-body interaction, is likely to intro- 
duce Triton on a large, highly eccentric orbit. Distant satellite orbits are perturbed primarily 
by the Sun, inducing precession of angular variables and oscillations in eccentricity and incli- 
nation, with minor semimajor axis variations. The two important periodic perturbations are 
those associated with 1/2 of the planet's orbital period ("evection") and 1/2 of the preces- 
sion period of the argument of pericenter u ("Kozai behavior"). For early Triton, evection 
was first noted by Goldreich et al. (1989); conservation of orbital angular momentum during 
tidal evoution implies initial pericenter of 7 Neptune radii (Rn] Triton's present orbit has 
a = 14Rn), and evection-induced oscillations in e produced minimum pericenters of 5Rn- 

However, Goldreich et al. (1989) ignored Kozai oscillations in e, which must have been 
present if Triton's inclination i at closest approach was same as now (i = 157°, measured 
with respect to Neptune's equator). Kozai osillations require that when u = 90° or 270°, 
both e and i = 157° are at a maximum (e and i oscillate in phase for retrograde satellites). 
Since almost all tidal dissipation occurs during this high-e phase of the Kozai cycle (when 
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the pericenter distance is smallest), this 157° inclination will be conserved as the maximum 
one for the Kozai cycle, while the minimum one (coinciding with u = 0° and 180°) will be 
affected by dissipation. Using more complete treatment of tides, Chyba et al. (1989) show 
that tides raised on Triton could not have led to a significant inclination change as long as 
it was in Cassini state 1. Trapping in the much less likely Cassini state 2 would have caused 
a rapid increase in its inclination (i.e., closer to 180°), rather than a decrease (Jankowski et 
al. 1989), so Triton's inclination relative to the local Laplace plane was always < 157°. This 
assumes orbital evolution slower than the nodal precession, preserving inclination relative to 
the local Laplace plane, which would initially be close to Neptune's orbital plane but would 
subsequently move closer to Neptune's equator. This condition is true for the tidal model 
(Goldreich et al. 1989) but not the gas-drag model (McKinnon & Leith 1995). Fig. 1 presents 
two illustrative short orbital evolutions. Using initial conditions of a = 2 x 10 7 km= 800-Ra? — 
0.18Rn,e = 0.982, u = 90°, and varying inclinations, we re-create two possible post-capture 
orbits for Triton. Both evection (P ~ 80 yrs) and Kozai oscillations (P ~ 1000 yrs) are 
clearly visible in the evolution of the inclined orbit, while the % ~ 180° case, shows only 
evection-related oscillations (whose amplitudes are in agreement with results from Goldreich 
et al. (1989) Fig. 2 for a = 800/iV). However, if tidal dissipation alone evolved Triton's 
orbit, only the inclined case can lead to present inclination of 157°. 

This conclusion points to a paradox. Goldreich et al. (1989) modelled the tidal orbital 
evolution with the standard relation: 



where a, t, q, n, k 2 and Q are respectively the semimajor axis, time, pericenter distance, mean 
motion, tidal Love number and tidal disspation factor, and Rt = 1353 km is Triton's radius. 
We numerically averaged Eq. 1 over several full Kozai periods, based on the output of the 
integration shown in Fig. 1. Using k 2 = 0.1 and Q = 100 for Triton (Goldreich et al. 1989), 
the resulting timescale for 1/e reduction in a is 3.5 Gyr for the inclined orbit and 0.11 Gyr 
for the near-coplanar orbit (cf. Goldreich et al. Fig.l). The requirement of an inclined orbit 
means that tides alone are not capable of circularizing Triton's orbit. 



Neptune's satellite system predating Triton's capture was likely similar to that of Uranus, 
as the planets are virtual twins in many important characteristics. Uranus possesses five size- 
able satellites at 5Ru < a < 23Ru, with total mass f«40% of Triton's. It is likely that Triton's 
capture led to destruction or ejection of any existing regular satellites outside 5R^. 



1 da 
a dt 
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(1) 



3. Interaction with Possible Regular Satellites 
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To explore the fate of these moons, we integrated five massless particles under the 
influence of Triton on the i = 157° orbit of Fig. 1. The innermost satellite was put at 
ai = 4.8Rn, and the others on Oj = 1.5aj_i; all the orbits were circular and in the plane 
of Neptune's equator (using current values of Neptune's obliquity and J 2 ). We found that 
the orbital crossing of the pairs 4-5 and 3-4 only a few centuries into the simulation. To 
calculate the collisonal timescale, we assigned our three outermost particles the radii of 
Umbriel, Titania and Oberon, and we took the satellites to be uniformly distibuted over 
all radii a(l — e) < r < a(l + e), and latitudes —i<(3< +i. Under these (admittedly 
rough) approximations, we estimated the collision probability for every output interval in our 
simulation and kept track of the cumulative collision probability for each pair. In both cases, 
the cumulative probability reaches unity at about 1000 years. While our approximations 
likely underestimate collision timescales, a more sophisticated method would not raise the 
result more than a factor of a few. Due to its long orbital period (about 7 yrs), most of 
it spent far from Neptune, eccentric Triton would be unlikely to collide with any of the 
satellites on such short timescales. 

The collision would likely destroy both moons and create a debris ring around Neptune. 
This debris would rapidly disrupt the remaining regular satellites, grinding them down to 
small pieces incapable of destroying Triton. The resulting massive disk out to about 20R N 
would interact strongly with Triton, whose perturbations would prevent re-accretion of the 
disk particles from taking place (cf. Banfield & Murray, 1992). The mutual interaction of 
the retrograde Triton and prograde disk, either through tidal torques or direct collisions, 
would cause decay of Triton's angular momentum; here we focus on the latter. 

To find the approximate evolution timescale of Triton's orbit due to the debris disk, 
we modeled the disk as having uniform surface density and extending out to 20R N , and 
then integrated Triton's orbit giving it a instantaneous "kick" every time it goes through 
Neptune's equatorial plane within 20R N . Like McKinnon & Leith (1995), we find that 
Triton's lifetime depends largely on its inclination: inclined orbits survive much longer than 
coplanar ones. Our timescales for the decay of orbits with % = 157° at to — 90° are on the 
order of a few times 10 5 yrs, while the orbit with % = 175° at u = 90° would have an e- 
folding time of a few times 10 4 yrs. While our disk model and resolution of the passages are 
crude, the robust conclusion is that the debris-drag timescale for orbital evolution is much 
shorter than the tidal one. Although our disk is much less massive than that of McKinnon 
& Leith (1995) (who use a minimum mass nebula) Triton can still evolve quickly due to 
deeper penetration into the disk. Finally, since the disk is less massive than Triton, there 
is a natural end to this stage of the evolution, so there is no need for specific mechanisms 
ensuring Triton's survival. 
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If Triton's evolution had been this fast (10 4 — 10 5 yrs), then the assumption of con- 
servation of inclination relative to local Laplace plane is incorrect, as the nodal precession 
timescale (at the boundary between an inner region dominated by Neptune's oblatneness 
and an outer region dominated by the Sun) is longer than 10 5 yrs (the transition distance 
depends on J 2 as the fifth root, making the influence of the additional material from the 
debris disk modest). Therefore, the pre-transition inclination of Triton could have been 
anything from 128° to 174° if the evolution was fast. Because passages through the disk 
would at that point (a ~ 80 — IOORn) happen for a wide range of values for u (and not only 
uj ~ 90° when e is high), the constraints on the initial inclination are very weak. 



4. Interaction with Irregular Satellites 

Nereid and the other five now-known irregular satellites of Neptune (Holman et al. 2004) 
can help us constrain Triton's post-capture orbit. During its early orbital history, Triton's 
orbit would intersect those of the irregulars, making gravitational scattering possible. There 
is reason to think that Neptune's irregulars had their orbits modified after capture: almost 
all irregular satellites of the other three giant planets have pericenters q around 100 and 
200 planetary radii, in case of direct and retrograde satellites, respectively. Four outer 
neptunians do not approach the planet closer than about 350Rn, while S/2002 Nl and 
Nereid come within about 50Rn', the latter two objects will be addressed in the next section. 

Can large-g irregulars result from Triton's passage? Fig. 2 shows six simulations using 
a symplectic integrator that includes Triton evolving under the influence of the disk while 
interacting with a swarm of test particles representing "original" irregulars. The 2044 test 
particles were originally put at a = 100 — 1820-Rjv, i = — 180°, u = 90° and with fixed 
q = 100 Rn (for direct) or q = 200 Rn (retrograde orbits). Triton's initial conditions were 
varied (see Fig. 2 caption). The curves in Fig. 2 give the fraction of surviving bodies over 
the 1 Myr integrations, normalized to the number of survivors when Triton is not present 
(some particles are inherently unstable). 

Fig. 2 shows that a significant number of irregulars survives Triton's passage only if 
the duration of the interaction is below ~ 10 5 yrs. In such cases, Triton's apocenter drops 
below the inner edge of the irregular population, and the depletion stops. In these models 
such rapid evolution occurs only for i max = 175°, for then Triton's pericenter is almost 
continuously within the disk. An orbit with this maximum inclination will occur for a large 
fraction of e > 0.9 post-capture orbits. Also, requirement for r < 10 5 yrs practically excludes 
tidal evolution independently from arguments in section 2. 
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Fig. 3 compares the semimajor axes and pericenters of Nereid and the five other irreg- 
ulars with the surviving test particles in the simulation represented by the topmost curve in 
Fig. 2. The scattered particles bridge the gap between initial conditions corresponding to 
other planetary irregulars and Neptune's. However, the mechanism leaves many satellites 
at q — 200 — 400i?jv. Interestingly, Nereid's apocenter distance seems to separate the bulk 
of the surviving particles and the observed moons. To test the importance of scattering 
by Nereid on the irregulars, we ran two simulations similar to the ones already described, 
except that now non-evolving Nereid was the scatterer and the simulations lasted 10 6 and 
10 7 years, respectively. There was a noticeable effect on the stability of particles, and some 
particles in the center of the "stable zone" were observed to escape. Nereid's escape speed 
(vesc = 0.25v or b at apocenter) apparently makes it an efficient scatterer, perhaps allowing the 
moons on orbits permanently out of Nereid's reach to dominate the surviving population. 



5. Early Dynamics of Nereid- Type Objects 

Nereid orbits between 55 and 385-Rjv, and has been previously suggested to be a regular 
satellite scattered to its current orbit by Triton. We report preliminary findings on the secular 
orbital evolution of Nereid-like objects just after the end of close encounters with Triton. By 
"Nereid-like" we mean that the precession is dominated by Triton, rather than the Sun. The 
radius of the zone within which the planetary oblateness dominates solar perturbations is 
r c = [2 J 2 (Mjv / Msun) R% a N] 1 ^ 5 > where Mn and Ms un are the masses of Neptune and the Sun, 
and a 7v is the planet's semimajor axis. We estimate the importance of Triton by putting 
J 2 = {M T R 2 N )/(2M N a\), where subscript T refers to Triton. If we take a T — l/2r c , we 
find r c = (0.25M T /M Sun ) 1/3 a N ~ 8 x 10 6 km. Therefore, "Nereid-type" objects are those 
orbiting within about 7% of Neptune's Hill sphere (for Nereid, a = 5.5 x 10 6 km= 220-Rat). 

The most important secular interaction between two non-intersecting orbits is through 
the quadrupole term in the interaction Hamiltonian (i.e. one containing a\/a\] subscripts 
1 and 2 refer to inner and outer body, respectively). The disturbing function for the outer 
body is (Innanen et al. 1997): 

R = ™^ Qr [2 + 3e 2 T - (3 + 12e 2 - 15e 2 cos 2 uo T ) sin 2 i 2 \. (2) 
8o 2 

where b 2 = a 2 (l — e 2 ,) is the outer body's semiminor axis, and and lot are those of 
Triton. If we assume Triton's orbit precesses much more slowly (due to Neptune's J 2 ) than 
the outer body's (due to Triton and the Sun), we can replace ujt with Q 2 , the longitude of 
the ascending node of the outer body, measured from Triton's pericenter, in Triton's orbital 
plane (the angles are geometrically identical). The expression above has no dependence on 
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oj 2 , so e = const. Using Lagrange's eqations, expressions for the evolution of i 2 and Vt 2 obey: 



where r = tMp/(b\n T )- Putting e T = 0.6 (which puts Triton close to intersection with 
Nereid, assuming current pericenters) , we numerically integrated Eqs. 3 and 4 for initial 
conditions of fi 2 = 90° and i 2 = - 180° (Fig. 4). Suprisingly, objects with % = 40° - 140° 
can have nodes librating around 90° (or 270°, as the results for Q 2 > 180° are symmetric), 
and their inclination (relative to Triton) can change from direct to retrograde. We tested this 
model using a numerical integration, in which Triton had its current % and q but ex = 0.6. 
Test particles were placed outside its apocenter, on e = orbits in Neptune's equator plane. 
Several particles show librating Q 2 and repeatedly change their sense of revolution (Fig. 4). 

Inclinations in Fig. 4 are relative to Triton's orbit, which has % = 157° relative to 
Neptune's equator, which itself has an obliquity of 29° relative to the planet's orbit. Objects 
in libration island can end up with a wide range of inclinations relative to the Sun. In 
particular, bodies starting within the equatorial plane of Neptune can be shuffled through 
the libration island to both direct and retrograde orbits. 

One recently discovered irregular, S/2002 Nl, has a = 1.6 x 10 6 km, q > 1.3 x 10 6 km, 
and i = 134° (Holman et al. 2004). Grav et al. (2004) report that Nl's colors are similar to 
that of Nereid and suggest that they are members of a collisional family. This is implausible 
if Nl was produced at the present epoch, but is consistent with both objects being pieces of 
circumplanetary debris ejected by Triton and mixed by this process. 



We propose that, after its initial capture by Neptune (via any process), Triton strongly 
pertubed all of the pre-existing satellites, both regular and irregular. Regular satellites 
rapidly collide with each other, creating a debris disk to be subsequently swept up swept up 
by Triton. The debris drag evolved Triton's orbit rapidly enough (< 10 5 yrs) to preserve 
some of the irregular satellites, scattering some to high q. Bodies in the inner part of the 
Hill sphere later suffered Triton's secular perturbations, which in some cases "flipped" their 
inclination between direct and retrograde. The closest irregular satellites were subsequently 
depleted through encounters with Nereid, producing the distribution observed today. 




(3) 



di 2 15 



e\ shu 2 sin(2fi 2 ), 



(4) 



6. Summary 
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Fig. 1. — Pericenter distances of Triton on large and eccentric post-capture orbits, with 
imin = 157° (at oo = 90°; solid line) and i = 180° (dotted line). When the pericenter 
is furthest from neptune, the inclination (with respect to the Neptune's orbit) is at the 
minimum value {i m in = 110° for i max = 157°). 
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Fig. 2. — Evolution of the fraction of surviving test particles over time. The labels on the 
graph give time Triton takes to evolve to half its initial a in each simulation. Dashed lines 
plot evolution with Triton's initial io = 157° and a = 10, 20, 30 x 10 6 km (right to left), while 
the solid lines plot simulations with i = 175° with a = 5, 10, 20 x 10 6 km (top to bottom). 
Neptune's hill sphere is 110 x 10 6 km in radius. 
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Fig. 3. — Semimajor axis-pericenter plot for simulated test particles (initial condi- 
tions:asterisks, final orbits:pluses) and real satellites (large boxes). Nereid's apocenter dis- 
tance is shown as a dashed line. 
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Fig. 4. — Evolution of early-Nereid-type orbits in i-Q space: dotted curves represent results 
of an analytic model while the solid curve shows the history of a numerically integrated 
particle [0,2 = 1.78ot; gt = 0.6 in both cases). The asterisk marks the initial condition. 
The inclination is measured with respect to Triton's orbital plane. The phase space for 
Q = 180° - 360° is identical, with another libration island at Q = 270°, % = 90°. 



